Indirect control of the asymptotic states of a quantum dynamical semigroup 



Raffaele RomancQ 
The Abdus Salam International Centre for Theoretical Physics 

Strada Costiera, 11 1-34014 Trieste, Italy and 
The Department of Theoretical Physics, University of Trieste, 
Strada Costiera, 11 1-34014 Trieste, Italy 

In the dynamics of open quantum systems, the interaction with the external environment usually 
leads to a contraction of the set of reachable states for the system as time increases, eventually 
shrinking to a single stationary point. In this contribution we describe to what extent it is possible 
to modify this asymptotic state by means of indirect control, that is by using an auxiliary system 
coupled to the target system in order to affect its dynamics, when there is a purely dissipative 
coupling between the two systems. We prove that, also in this restrictive case, it is possible to 
modify the asymptotic state of the relevant system, give necessary conditions for that and provide 
physical examples. Therefore, in indirect control schemes, the environmental action has not only 
a negative impact on the dynamics of a system, it is rather possible to make use of it for control 
Oh purposes. 
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Introduction 



The study of quantum mechanical systems is relevant not only for a deep understanding of the fundamental physical 
laws, but also for its potential applications In particular, it has been proved that the use of quantum-based 
technologies would highly increase the performance of a computational device ■ This is due to the mathematical 
' structure of the theory, leading to peculiar features for the microscopic systems, absent in the macroscopic, classical 
world. Among them, the most relevant (and not completely understood) is represented by the quantum correlations 
CX' known as entanglement, typical of quantum systems, whose complete characterization has been given only in the low 
dimensional cases. These correlations are the key of many recently proposed protocols, as teleportation ^3'] or quantum 
7-H . cryptography 

' However, quantum systems are fragile: their relevant features are usually degraded by their interaction with the 
, external environment, leading to irreversibility, dissipation and decoherence. Therefore closed systems, described by 
the Schrodinger equation, are only approximations of real systems, that are necessarily open since they exchange energy 
and information with the external world. To account for this, the standard dynamical model is given by a quantum 
dynamical semigroup, that is a one-parameter family of Markovian (i.e., satisfying the semigroup property) completely 
positive maps, transforming pure states into mixtures and destroying quantum coherence [5|, l6(- Irreversibility is due 
to the fact that quantum dynamical semigroups are contractions in the state space of the system: this highly reduces 
the ability of manipulation on the system, since many transitions become forbidden. The extremal case is represented 
by the so-called uniquely relaxing semigroups, where there is a unique asymptotic state for every initial state of the 
system. Therefore, in order to fully make use of the potentialities of quantum systems, it is fundamental to understand 
the mechanisms leading to dissipation and to conceive methods to counteract them, and, more in general, to control 
■ the dynamics of the system. 

With these motivations in mind, several solutions to the problem have been proposed. The study of the symmetries 
of the system-environment coupling has lead to the notion of decoherence- free subspaces or subsystems, that are 
uiet places, unaffected by decoherence, where to encode the relevant information (for a review on the topic, see 
). From a more active perspective, a theory of quantum control has developed, dealing with the effect of external 
manipulations that can be performed on the system. In the standard setting [coherent control), the control parameters 
enter the Hamiltonian of the system, for example through external fields c oup led to t he sy stem (for a geometric control 
perspective, in line with that developed in this work, refer to [1, [13, El, [H, [H, Although this is the most 

natural way to introduce external actions in the dynamics of a system, its ability to fight the unwished effects of the 
irreversible dynamics is limited [3], unless some information about the state of the system is collected, and then used 
to update in real time the controls. This information is usually obtained via an indirect continuous measurement 
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and, because of this, the master equation describing the system becomes stochastic. This quantum feedback scheme 
represents a promising approach in many interesting cases [HI, [l(| [l7| • 

A different approach to quantum control has been recently discussed, in which the control does not enter through 
the Hamiltonian of the system, but it is rather obtained by means of an auxiliary system (ancilla) [IH, QjJUil) that 
can be manipulated and put in interaction with the relevant system, and finally discarded at the end of the procedure. 
The ability of driving the target system through the auxiliary one is determined by the correlations between them. 
The interest in this indirect control method is twofold. It is complementary to the coherent control approach, that is, 
it can be applied to experimental setups where the coherent control technique is not appropriate. Moreover, in the 
indirect control approach the environment does not only represent a source of noise, it can also be used for control 
purposes. In fact, it has been proved that the environment can correlate two systems immersed in it via a noisy 
mechanism, not only destroy their existing correlations [2l|, [22j • This mechanism can be used to obtain total control 
of the target system even if there is not a direct interaction between the two parties [23| . 

In this contribution we prove that, in the indirect control scheme, the environment induced correlations can be used 
to manipulate the asymptotic states of the target system. This result gives further evidence that, in this framework, 
it is possible to engineer the environmental action to get controllability. The work is organized as follows: in Section U] 
we summarize some standard results about quantum dynamical semigroups and their stationary states, in Section HT1 
we provide necessary conditions for the indirect manipulation of the asymptotic states of the relevant system under a 
purely dissipative dynamics, in Section HTT1 we evaluate the stationary states and describe a concrete physical example 
of application, and finally we conclude in Section HVl 



I. QUANTUM DYNAMICAL SEMIGROUPS AND THEIR STATIONARY STATES 

In many situations (usually, when there is a weak interaction with the surrounding environment) it is possible 
to approximate the reduced dynamics of an open system S using a Markovian one-parameter family of completely 
positive maps {7*; t ^ 0}, satisfying the semigroup property 7*+ s = 7* ° 7s, with t, s ^ 0, with 

p s (i)= 7t [p s (0)], (1) 

where the state of the system S is given by the Hermitian, positive, unit trace operator p s (statistical operator), 
acting on the n-dimensional Hilbert space associated to S. Complete positivity is necessary in order to guarantee a 
physically consistent interpretation of the formalism when dealing with composite, entangled systems. The generator 
L of the dynamics is defined by p s = L[p s ], and it is possible to prove that it has the general form (the so-called 
Lindblad-Kossakowski form) 

L[ Ps ] = -i[H s , p s ] + ]T dj (f iPs f} - \{FjF u Ps }) , (2) 

where H 8 = H\ is the Hamiltonian of S, and the set {Fi;i — 1, ... ,n 2 — 1}, along with the n-dimensional identity 
operator, form a basis for the operators acting on the Hilbert space associated to S, satisfying TrFi = 0, and 
Tr(FiF^) = 5ij. The (n 2 - 1) x (n 2 - 1) matrix C = [c tJ ] (Kossakowski matrix) fulfills C t = C and C ^ 0, necessary 
and sufficient condition for the complete positivity of the dynamics [H, Q . 

In the following, S is a bipartite system, S = T + A, where T is the target system (to be manipulated) and A the 
ancilla. We assume that T and A are two copies of the same two-level system, separately interacting with the same 
environment, assumed to be spatially homogeneous, according to the Markovian dynamics ([2|). We further assume 
that H s = 0, since we want to study a purely dissipative dynamics. To model this system, it is sufficient to consider 
the basis {Fi; i = 1, . . . , 6}, given by the local operators Fi = <?% ® I for i = 1, 2, 3 and Fi — I<E> 0^-3 for i = 4, 5, 6, 
where cr^, i = 1,2,3 are the Pauli matrices. We consider the standard representation of these operators in which 173 
is diagonal. The 6x6 matrix C has the form 



C 



A B 
St A 



(3) 



where A — A^ is the Kossakowski matrix for the system T (or A) alone, and B represents the dissipative coupling 
between the two parties. A and B are 3x3 blocks. The form is not the most general joint Kossakowski matrix, 
as we have assumed that the two parties interact separately with the environment, and that the two local dissipative 
contributions are equal (homogenous environment). We will limit our attention to models well described by ([3]); 
moreover, for simplicity, we will further assume B = B' . This choice produces a significative simplification in the 
treatment and it is still of great phenomenological interest. 
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The first term in the right hand side of ([2]) represents the coherent part of the evolution and it generates a group of 
reversible, unitary transformations. The second term generates the irreversible dynamics, according to the matrix C 
whose entries depend on the microscopical details of the interaction between system and environment. It also leads 
to the appearance of attractors in the state space of S, and consequently relaxation to equilibrium of the states of 
the open system. A stationary state for the dynamics, pf, is defined by the condition on the generator = 0. 

Since the dynamics is linear in the state p s , it is possible to fully characterize the asymptotic fate of the system 
by studying the eigenvalues of the dynamical matrix appearing in the coherence vector representation of @ [HI- 
Although this treatment is very general, it is not suitable for the purposes of this work. We will rather refer to some 
necessary conditions for the existence of stationary states, and for the convergence of p s (t) to them, given in terms of 
the operators {Vi \ i} appearing in the diagonal form of (O, 

Hps] = E(^^ - \i v ^Ps})- (4) 

i 

These conditions are expressed by the following theorem [25[, that has been adapted to the present context. 

Theorem 1 Given the quantum dynamical semigroup Q), assume that it admits a stationary state po of maximal 
rank. Defining At = {Vi,V^ ;i}' , the commutant of the Hamiltonian plus the dissipative generators, the following 
conditions hold true: 

1. If At = span(I), then po is the unique stationary state. Moreover, if {Vi;i} is a self-adjoint set with {Vi',i} = 
span(I), then for every initial condition p s (0) 

lim p s (t) = p . 

t — >+oo 

2. If At ^ span(I), then there exist a complete family {P n ;n} of pairwise orthogonal projectors such that Z ~ 
Ai D At = {P n ;n} . If {Vi\i} — Ai, two extreme cases together with their linear superpositions may occur. 

i. If Z = At, then for every initial condition p s (0) 

lim p.(t) =^Tr(P nPs (0)P») T f^° P " ■ 



ii. If Z — At , then for every p s (0) 



lim ft (t) = VP nft (0)P„. 

t— »+oo * — * 



Therefore, the stationary states of a quantum dynamical semigroup can be characterized by means of the algebras 
Ai, At' , and Z, if a maximal rank stationary state po is available. These quantities are evaluated in the next section, 
depending on the form of the matrix C . 



II. RELEVANT ALGEBRAIC QUANTITIES 



Following Theorem [I] we need to write C in diagonal form in order to find the operators Vi appearing in (|U). This 
is achieved by means of the unitary transformation U such that 



UCU^ = diag(A t ,i = 1,. ..,6), 
where \ are the eigenvalues of C. U has the form 



U 



V2 



u u 
-u ii 



(5) 



(6) 



and U, U are unitary transformations such that 

U(A + B)W = diag(A+,i = 1,2,3), 
U(A-B)fr = diag(Ar,z = 1,2,3). 



(7) 
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The eigenvalues of C are ordered as A; = for i = 1,2,3 and Xi — X i _ 3 for i — 4,5,6. Comparing the generator 
forms © and (d]), and using the notation U = [iiy], we have 

6 

k=l 

Following ([6]), it is possible to write 

. fl®^ + ff,®l, i= 1, 2, 3 

Vi={ (9) 
I (g) (7j_3 - (Tj_3 ® I, i = 4, 5, 6 



where we have defined 



3 3 



fc=l fe=l 

and we used the notation U = [uy], J7 = The operators in (|10p satisfy Tro^ = Tro^ = and Tr(cT'icrt) = 

Tr(CTi(rj) = &y. They are self-adjoint if and only if the unitary operators U and U are orthogonal. 
The commutant of Theorem [T] can be expressed as 

M = {V i ,Vt;i\\ i ^0}' = f| {VS,^}', (11) 

where only non-vanishing eigenvalues Xi have to be considered, otherwise the corresponding Vi do not appear in the 
generator ([4]). Moreover, for a given i, 

{Vi,V?}' = {v\ve{Vl}',vie{Vi}'}, (12) 

therefore we can limit our attention to the sets {Vi} . We find convenient to consider separately the two kinds of 
contributions defined in ©. To begin with, we consider a fixed index i such that 7^ 0, and assume that the 
corresponding ex, is non-singular. In this case it can be written as 

Gi = fxiRidzR' 1 (13) 

where it is possible to choose Ri = R^ \ and 

m? = Ek-) 2 - ( l4 ) 

3 

Since 1® &i + &i®I = jlilZi (I ® 0-3 + 03 ® T)7Zi , with TZi — Ri® Ri, it follows that 

{I ® CTj + ct 4 (g) I}' = 7^{I cr 3 + cr 3 O 1}'^ (15) 
and then, after the explicit computation, 

{V}' =span(I®I,l(g)CT i ,a- i <g)I,CT i ®a i) fi + ,Ar), (16) 
having defined the additional operators 

fl + = <T\ ® <Ti + (T 2 ® CT 2 + CT3 ® CT3 , 

A," = ^ 4 (CTi®Cr 2 -CT2«)Cri)7e. i . (17) 

Notice that, in general, the operators in the right hand side of (fT6|) are not self-adjoint, nor orthogonal each other 
in the Hilbert-Schmidt metric, since the transformation TZi is not unitary (equivalently, self-adjoint). However, if the 
coefficients u*,-, j — 1,2,3, are real, (7j is self-adjoint and TZi unitary. Consequently, in this case the basis of {Vi} is 
made of Hermitian, orthogonal operators. 
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The commutants {Vi} are completely characterized for i = 1,2,3. Finally, {Vi, V^} can be found by considering 

{Vi,V?}=\ (18) 
[ span(I ® I, il + ), otherwise. 

The corresponding sets for i = 4, 5, 6 can be found by applying the same procedure to <7j, assuming that A^~ 7^ 0. 
The result is 



where 



and Si = & ® & , with 



where Si = , and 



Finally, in this case 



{V;} = span(I ®I, I ®3i, 3* 81,^®^, fir, A+), (19) 



SI,, = 5j((Ti ® (Ti - cr 2 ® cr 2 )iSi, 

A+ = Si(ai ®<j 2 +a 2 ® o-i (20) 



= fi i S i a 3 S i \ (21) 
A? = EW- (22) 



{^V- 1 } = ^ (23) 
[ span(I®I), otherwise. 

If cf-j (or (Ti) is singular, the previous computations are not longer valid. In this case, the commutants must be 
evaluated by direct computation and it is not possible, in general, to express their structure in a compact form. 

We have all the ingredients to evaluate the contribution related to the dissipative generators Vi in (fTTj) . We find 
convenient to denote by n + and n_ the number of non-vanishing eigenvalues of the type A + and A - respectively. 
The non-trivial cases are summarized below, with the corresponding relevant algebras and set of projectors {Pf,i}, 
described in theorem [1] to be used to construct the set of stationary states. For further reference, the projectors Ilfc, 
k e { — , +, 1, . . . ,4}, are defined as 

n fc = [7r£], tt^ = 5 ik S jk , k = 1, . . .4; 

n_ = -(i<g>i-fi+), n + = i®i-n_. (24) 

A. Case 1 

It is characterized by n+ = 1, A = A T , B = A. We notice that A = B is equivalent to n_ = 0. The commutant 
is given by M = span(I ® I, I ® o-j, <Xj ® I, (Xj ® 5j, fi + , A~), and Z = .M' = span(I ® 1, cjj ® crj,! ® cr% + <Ji ® I). The 
projectors are given by 

Pi = T^iirfti, p 2 = n t u 4 n u 

p 3 = ^(n 2 + n 3 )^. (25) 

B. Case 2 

It is characterized by n + — 1, A ^ A T , £? = A, or rather n+ > 1, B = A. In this case Z = M. = span(I ® I, fi + ), 
and there are only two projectors: 

Pi = n_, p 2 = n+. (26) 
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C. Case 3 

It is characterized by 714. = n_ = 1, A = A T , B = aA, a € K \ { — 1, +1}< We observe that [j4, B] = 0, thus it is 
possible to choose U = U. Moreover, B = aA implies = <7£ for the index £ such that A^ 7^ and A^~ 7^ 0. Finally, 
Z = A4 = span(I ® I, ct, ® it,,! (8 o"j, <7i (g) I), and the projectors are given by 

/'., K.U/R... ./ 1 1. (27) 

In all the remaining cases M. = span(I ® I), part 1 of Theorem [T] applies and the maximal rank stationary state is 
unique (if there is one). Therefore, the aforementioned cases are necessary conditions for the indirect manipulation 
of the asymptotic state of the target system T via the auxiliary system A. 



III. STATIONARY STATES 



We separately explore the non-trivial cases described in Section [Til Following Theorem [TJ if a stationary state po 
whose eigenvalues are all non- vanishing can be found, it is possible to build the whole family of stationary states , 
by using the projectors {Pi; i}. Finally, the corresponding stationary state of the target subsystem can be obtained 
from 

p?=Tr A p?, (28) 

that is by a partial trace over the degrees of freedom of the auxiliary system. We consider two different choices for the 
initial state p s (0), depending on wether there are initial correlations or not. As the first choice we take into account 
the product state 

p s (0) =pr(0)®p A (0), (29) 

where Pt(0) and pa(0) are arbitrary states for the two subsystems, that will be written using a Bloch vector repre- 
sentation as 

1 3 

pT(0) = -(l + 5>£cx fc Y (30) 



fe=i 



with real coefficients p\, and analogously for pa(0), with real coefficients pf. The choice (|30"|) refers to initially 
uncorrelated systems, that will in general couple during their joint, dissipative evolution, because of the off-diagonal 
block B in the Kossakowski matrix C. Alternatively, we consider the pure initial state 



p.(0) = M(V|, \4>) = vP\ T) ® I T) + ^~p\ i) ® I I), (31) 

where Pel, and | f), I -L) are the +1, respectively —1 eigenvalues of the operator 03. This state is entangled if 
P ^ 0,1, and it is maximally entangled if P = |. It is not an arbitrary entangled state, nevertheless it can be used 
to test the impact of an initial quantum correlation on the manipulation of the stationary state of T. 

Although their algebraic structures are different, Cases 1 and 3 lead to the same results. Since Z = A4', in Case 1 
there is not need of po, whereas the simplest maximal rank stationary state in Case 3 is given by the maximally mixed 
state po — I <8 I. If there is not correlation in the initial state, it is not possible to manipulate the stationary state 
of the system T by means of the ancilla A. In fact, the coefficients of Bloch vector representation of p^ , denoted by 
pf 3 , i = 1,2, 3, depends only on pr(0): 

PT = u ii (Pi^ei ~ P2 u& + P3 U ^) 

P? = -u^2[pTu^- pIu^ 2 + P^Utf) (32) 
Pf 



where £ <E {1, 2, 3} is such that A^ 7^ in Case 1, \f 5^ and A^ 7^ in Case 3. If we consider the (possibly entangled) 
initial state (|31[) . the dependence on P is apparent: 

p? = (2P-l) HlH3 

p? = -(2P - l)u i2 u i3 (33) 



p? = (2P-l)u 2 
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Therefore, at different correlated initial states there correspond different stationary states pjt. Manipulations of the 
asymptotic states of the target system are possible only when there is an initial correlation between T and A. 

In Case 2, since the expression of the stationary state p 1 ^ is more involved, we prefer to present a concrete example 
in which both uncorrelated and correlated initial states allow indirect manipulations of the asymptotic states. A 
simple case is given by the choice 



A = B = 



a ib 
-ib a 
0a 



(34) 



with the condition a 2 — b 2 expressing the complete positivity of the evolution. In this case, the maximal rank 
stationary state is found to be 



PO = 7(l(g>I+-(I<g)0-3 + cr3<gSl)+ (-) <7 3 <E> a 3 ) 
4 V a Va/ / 



6\2 



4 \ a \a/ 

and the asymptotic state of T for the uncorrelated initial state has components 



(35) 



pT = o 
p? = o 



(36) 



P? = 7-1 3- 



ET A 
Pk Pk 



k=l 



where 



For the correlated initial states we get 



ab 



3a 2 + 26 2 ' 



(37) 



PT 
Pi 



= 
= 



(38) 



pf = 4r 1 + V /P(1-P) 



Therefore, in both cases it is possible to manipulate p^- We stress that it is possible to drive the asymptotic state 
of T through the initial state of A even if the two system are initially uncorrelated, and there is not a Hamiltonian 
coupling between them. 



IV. DISCUSSION AND CONCLUSIONS 



We have explored the asymptotic performance of the indirect control method when both target and auxiliary 
systems are two-level systems, and they evolve under a purely dissipative dynamics. We have assumed that the two 
systems interact separately with an homogeneous environment, leading to a particular form of the Kossakowski matrix 
C for the composite system. We have found that the conditions expressed in Case 2 are necessary conditions for the 
indirect manipulation of the stationary state of T through the initial state of A when the initial state is a product 
state. We have given a numerical example in which this dependence is explicit. We have also proved that, in all the 
non-trivial cases considered, an initial entanglement between the two systems is effective for control purposes. 

Initial states with a different correlation between the two parties produce different stationary states for a reduced 
subsystem whenever some correlation survives to the decohering action of the environment. This is also true for more 
general models than the one described in this contribution. 

For initially uncorrelated states, the dissipative evolution has to provide the necessary entanglement, that has to 
be preserved in the large time limit (for the asymptotic entanglement in a quantum dynamical semigroup with purely 
dissipative evolution, see the results presented in [26j]). Therefore, the ability of varying the stationary state of T 
is a controlled dissipative mechanism. In this sense, in the indirect control approach the environmental action can 
be considered as a resource. This kind of behavior has already been observed when dealing with accessibility and 
controllability of a pair of qubits immersed in a bath of decoupled harmonic oscillators, in an exactly solvable model 
[23| . Therefore, it is not an artifact of the Markovian nature of the evolution. 
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Eq. (|37|) gives a limited ability of manipulation of p^?. However, the case here discussed is intended to represent 
an example of explicit dependence, not a complete treatment of the asymptotic reachable set. Moreover, the relevant 
case A = B is important in concrete experimental situations, for example in the study of the resonance fluorescence 
[27L [28| , or in the analysis of the weak coupling of two atoms to an external quantum field [29[ . 
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